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The topological susceptibility in 'full' (UK)QCD. 
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We report first calculations of the topological susceptibility measured using the field theoretic method on SU(3) 
gauge configurations produced by the UKQCD collaboration with two flavours of dynamical, improved, Wilson 
fermions. Using three ensembles with matched lattice spacing but differing sea quark mass we find that hybrid 
Monte Carlo simulation appears to explore the topological sectors efficiently, and a topological susceptibility 
consistent with increasing linearly with the quark mass. 



The latest UKQCD runs generate SU(3) field 
configurations using a Wilson gauge action cou- 
pled to n/ = 2 flavours of Wilson sea quarks, 
non-perturbatively improved such that the lead- 
ing order discretisation errors in spectral quan- 
tities should be 0(a 2 ). In dynamical simula- 
tions the lattice spacing, a, is influenced by the 
gauge coupling, (3, and the quark mass parameter, 
k. The ensembles studied represent three points 
on a trajectory in the (j3, k) space of (approxi- 
mately) constant a and physical volume (Ua 4 = 
16 3 .32a 4 ) §. [The lattice spacing has been de- 
fined from the static inter-quark potential, using 
ro = 0.49 fm 0].] Along such trajectories, how- 
ever, there is a change in chiral behaviour as seen 
in the pseudoscalar to vector mesonic mass ratio 
in Table |l]. [Dimensionless lattice quantities are 
denoted with circumflexes throughout.] 

We measure the topological charge density on 
the lattice using the field theoretic operator 

Q{n) = ~ x ^e% aT Tr U^(n)U aT {n), (1) 

symmetrised such that £ _ T = -£ +T etc. 
to improve the signal on moderately cooled lat- 
tices. The topological charge, Q, and susceptibil- 
ity, x, are then 
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Although formally in the weak coupling limit 
Q is related to that of the continuum by Q = 



a 4 Q + C(a 6 ), the prefactors of the higher order 
terms typically have momenta on the scale of 1 /a 
such that they are 0(1), and Q becomes increas- 
ingly dominated by ultraviolet noise near the con- 
tinuum limit. It is also subject to a multiplicative 
renormalisation that reduces the signal at the (3 
couplings accessible to simulation. 

An established, and well understood, proce- 
dure to extract the topological signal is to cool 
the configurations prior to measuring Q. This lo- 
cally smoothens the lattice fields to remove the 
ultraviolet fluctuations, and drives the renormal- 
isation constant to unity. We move through the 
lattice links in a staggered fashion, updating each 
Cabibbo-Marinari subgroup in turn so as to ex- 
actly minimise the SU(2) Wilson plaquette ac- 
tion. This action is the most local and thus is 
expected to do least damage to correlation func- 
tions on physical length scales. Updating every 
link once corresponds to one cooling 'sweep'. 

Such a cooling action, however, also destroys 
topological features when applied in extremis. 
Instanton-anti-instanton pairs in the vacuum are 
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5.20 

0.1350 

150 

4.58 (8) 
0.69 (1) 



5.26 

0.1345 

100 

4.58 (6) 
0.79 (1) 



5.29 

0.1340 

100 

4.45 (6) 
0.83 (1) 



Table 1 

Lattice spacing and pion to rho mass ratios 
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Figure 1. Normalised (Q(n c ) ■ Q(25)) versus the 
number of cools, n c , at (3 — 5.20. 

not stable minima of the action (even in the con- 
tinuum) and under cooling there is an attractive 
force leading to annihilation. Whilst this is an 
issue in measurements of instanton size distribu- 
tions, there in no net change in the topological 
charge and the susceptibility is thus stable. 

The lattice regularisation breaks scale invari- 
ance, and the Wilson action of an instanton is an 
increasing function of the core size, p. (Isolated) 
instantons shrink and disappear under cooling, 
leading to a net change in Q and \. Whilst such 
events may be monitored in relatively smooth 
configurations, it is still desirable to perform as 
few cooling steps as possible to expose the signal 
before making measurements. In Fig. ^ we plot 
the normalised correlation between the topologi- 
cal charge after a given number of cooling sweeps, 
n c , relative to that after n c = 25: 

(Q(n c ) ■ Q(25)) 
i((QK)) 2 + (Q(25))2) 

We find remarkable stability in Q from 5 cooling 
sweeps out to at least 25 cooling steps. 

One technical point of interest is the rate at 
which configurations become topologically inde- 
pendent under Monte Carlo updating; the topo- 
logical charge is related to the small eigenvalues of 
the fermion matrix and should be one of the slow- 
est modes to decorrelate. We show a Monte Carlo 
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Figure 2. Q(n c — 10) vs. simulation time in units 
of 40 HMC trajectories at f3 = 5.20. 

time history plot for our most chiral ensemble in 
Fig. ||, and we estimate the integrated autocorre- 
lation times in units of the 40 HMC trajectories 
between configurations in Table [|, albeit using a 
small ensemble, indicating pretty good decorrela- 
tion at these quark masses. We note that auto- 
correlation times are longer for the larger j3 en- 
sembles despite these being the less chiral. 

Also of issue is the ergodicity of the MC up- 
date. Measuring Q over the whole configuration 
we find the ensemble averages to be within ~ 1.5 
(relatively large) standard deviations of zero in 
Table and Fig. || shows a relatively Gaussian 
sampling of the topological sectors. 

In the continuum, the topological charge of a 
configuration is integer. On the lattice this is not 
so, nor is the charge measured using our operator 
on a configuration particularly close to integer. 
This can be attributed to the presence of nar- 
row instantons whose charge is significantly less 
than unity. Attempts can be made to calibrate 
and correct for this but we defer application 
of such procedures until || . We have both main- 
tained the charge as a non-integer for the pur- 
poses of calculating the topological susceptibility, 
and also rounded the value on each configuration 
to the nearest integer, and found that the topolog- 
ical susceptibilities obtained are consistent within 
statistical errors. 

In Table § we show our estimates for the topo- 




Figure 3. Histogram of Q lnt at n c = 10 on the 
/3 = 5.20 ensemble. 

logical susceptibility measured after 25 cooling 
sweeps. In the chirally broken, confining phase 
at low temperatures, the sea quarks induce an 
attractive interaction which leads to instanton- 
anti-instanton pairing and a suppression of the 
topological susceptibility as the dynamical quark 
mass is lowered: 
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where (f/'V') is summed over light quark flavours, 
and should be evaluated in the m q — > limit. The 
quark mass is not known a priori but may be re- 
expressed in terms of the pseudoscalar decay con- 
stant using the PCAC relation rriqfyip) = f%rn% 
such that for sufficiently chiral sea quarks, the 
topological susceptibility should be quadratic in 
the pseudoscalar mass and decay constant: 
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5.20 


0.71 (41) 


0.03 (56) 


7.5 (1.4) 


5.26 


0.67 (28) 


-0.47 (28) 


7.2 (1.5) 


5.29 


1.5 (1.2) 


1.10 (97) 


14.8 (2.3) 



Table 2 

Mean Q, autocorrelation time and the topological 
susceptibility at n c — 25. 
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Figure 4. ro(x mt ) 4 vs. \fr^m^ with a linear fit 



In Fig. ^ we plot the susceptibility versus the 
pseudoscalar mass in units of ro, and find the 
data are consistent with such a leading order fit 
through the origin. We may use the fitted slope 
of this graph to provide an estimate of the pseu- 
doscalar decay constant, finding it to be very low 
compared to the experimental value. 

Although such an estimate must be regarded as 
preliminary given the volume of data so far anal- 
ysed, it is clear from Fig. || that the heaviest of 
our sea quark masses yields a topological suscep- 
tibility that is in statistical agreement with the 
quenched value, and attempting to fit the lead- 
ing order chiral behaviour to this point is likely 
to lead to an underestimate of /„■; 0(rn^) terms 
are likely to be large, including the effects of not 
extrapolating (vpip) to m q — > 0. Ongoing analysis 
of further configurations at these and lighter sea 
quark masses should clarify the situation. 
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